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ABSTRACT 

Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades 
the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming 
region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding H ii region on 
the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration 
following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line- 
of-sight component of the magnetic field of about 3 pG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 
12 cm“*. The dust shell that surrounds the Rosette H ii region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal 
significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky 
orientation of the magnetic field in the Rosette’s parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. 
The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for 
a given orientation of the field. This comparison leads to an upper limit of about 45° on the angle between the line of sight and the magnetic field 
in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4 %. This field direction can reproduce the RM values detected in the 
ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5-9 pG. The present analytical model is able to 
reproduce the RM distribution across the ionized nebula, as well as the mean dust polarization properties of the swept-up shell, and can be directly 
applied to other similar objects. 

Key words. Polarization - ISM: magnetic fields - radiation mechanisms: general - radio continuum: ISM - submillimeter: ISM 


1. Introduction 

Stellar bubbles are formed from the combined action of ioniza¬ 
tion po'wer and stellar 'winds from OB stars. They are a major 
source of turbulent energy injection into the interstellar medium 
(ISM), sweeping up the surrounding gas and dust and modify¬ 

* Corresponding author: M. 1. R. Alves marta.alves@irap.omp. 

eu 


ing the magnetic fields (e.g. Ferriere 2001). Bubbles are seen 
as agents of triggered star formation, as the swept-up material 
may become unstable and fragment. It is therefore important to 
understand the influence of star formation and feedback on the 
structure of the magnetic held, and how the ionized bubbles are 
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shaped by the field. The all-sky polarization data from Planck^ 
open up new opportunities for studying the magnetic field struc¬ 
ture of such objects and their connection with the large-scale 
Galactic field. 

The interplay between the expansion of an ionized neb¬ 
ula and the action of the magnetic fields has been the subject 
of several studies, both observational (e.g. Pavel & Clemens 
2012, Santos et al. 2012, Santos et al. 2014) and numerical (e.g. 
Bernstein & Kulsrud 1965, Giuliani 1982, Ferriere et al. 1991, 
Krumholz et al. 2007, Peters et al. 2011,Arthur et al. 2011). 
These observational studies rely partly on interstellar dust grains 
as tracers of the magnetic field. Non-spherical grains spin around 
their axes of maximal inertia, which are precessing around the 
magnetic field lines. The grains emit preferentially along their 
long axes, thus giving rise to an electric vector perpendicular to 
the magnetic field. At the same time, dust grains polarize the 
light from background stars and since extinction is higher along 
their longest axes, the transmitted electric vector is parallel to 
the magnetic field in the plane of the sky (see e.g. Martin 2007). 
Observations of dust emission or absorption allow us to retrieve 
the plane-of-the-sky orientation of the magnetic field that per¬ 
vades interstellar matter. Results from high-resolution studies of 
H II regions have revealed regions of well-ordered magnetic field 
along the edges of the nebulae and magnetic field strengths of 
tens to hundreds of yuG. Santos et al. (2014) studied the star¬ 
forming region Sh2-29 in optical and near-infrared polarime- 
try and derived a field strength of around 400 //G. The authors 
used the Chandrasekhar-Fermi method (Chandrasekhar & Fermi 
1953), which relates the dispersion in magnetic field orienta¬ 
tion with turbulent motions of the gas, under the assumption of 
equipartition between magnetic and turbulent/thermal pressures. 
Such high values of the field strength reflect the ordered struc¬ 
ture of the magnetic field at sub-parsec scales, compressed by the 
expanding H ii region. The aforementioned dust polarization ob¬ 
servations of ionized nebulae, both in emission and extinction, 
have covered small regions around the objects, generally tens 
of arcminutes. Planck data allow us to study the magnetic fields 
probed by polarized dust emission towards large H ii regions em¬ 
bedded in their parent molecular clouds, and the diffuse medium 
surrounding them. 

Measurements of the magnetic fields towards H ii regions can 
also be performed through Faraday rotation of linearly polar¬ 
ized background sources, both extragalactic and Galactic (Heiles 
& Chu 1980). The plasma in the Hii regions rotates the plane 
of polarization of the background radio wave by an angle that 
is proportional to the square of the observing wavelength. The 
quantity derived from Faraday rotation observations, the rotation 
measure (RM), is directly related to the line-of-sight component 
of the magnetic field By weighted by the electron density of the 
ionized gas. Harvey-Smith et al. (201 1) studied the line-of-sight 
magnetic field towards five large diameter Galactic H ii regions, 
tens of parsecs wide, using RM combined with Ha data to esti¬ 
mate the electron density. They found By values of 2-6 /iG, con¬ 
sistent with those measured in the diffuse ISM through Zeeman 
splitting observations (Heiles & Crutcher 2005; Crutcher et al. 
2010; Crutcher 2012). This indicates that the RM enhancement 
observed towards H ii regions may be the consequence of a lo- 


* Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided hy two sci¬ 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci¬ 
entific consortium led and funded by Denmark. 
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Fig. 1. The Rosette Nebula, its molecular cloud Mon OB2, and 
their surrounding medium as seen by Planck at 353 GHz. The 
dotted lines give the approximate outer radii of the Rosette 
H ii region and the Mon Loop supernova remnant (described in 
Sect. 2). The northern Monoceros cloud Mon OBI is above the 
Galactic plane; part of the dust ring of A Ori is visible at the bot¬ 
tom right corner of the map; the northern end of Orion B is seen 
at (/, b) = (204°, -12°) (Dame et al. 2001). 

cal increase of the electron density. A similar study performed 
on the Rosette nebula by Savage et al. (2013) attributed the high 
values of RM to an increase in the field strength, as will be dis¬ 
cussed in Sect. 2. 

The goal of this work is to study the structure of the mag¬ 
netic field in and around an ionized bubble created by young 
stars. For this purpose we chose the Rosette nebula, adjacent 
to the Monoceros (Mon) OB2 cloud, for its close to spherical 
shape and for its relatively large size on the sky, 1°4, relative to 
the Planck 353 GHz beam of 4.'9. This paper is organized as fol¬ 
lows. We start in Sect. 2 by describing the Rosette nebula and 
its main features relevant for our study. In Sect. 3 we introduce 
the Planck and ancillary data used in this study, which are then 
analysed and discussed in Sect. 4. The interpretation of the ra¬ 
dio and submillimetre polarization observations is presented in 
Sect. 5, in light of a 2D analytical model of the magnetic field 
in a spherical bubble-shell structure. The main results are sum¬ 
marized in Sect. 6. The detailed derivation of our magnetic field 
model is given in Appendix A. 

2. The Rosette Nebula and the Mon OB2 molecular 
cloud 

The Rosette bubble is located near the anti-centre of the Galactic 
disk, centred on (/, b) - (206?3, -2?1), and is part of a larger re¬ 
gion called the northern Monoceros region. It is at (1.6 + 0.1) kpc 
from the Sun and its age is estimated to be about (3 + I)Myr 
(Roman-Ziiniga & Lada 2008). Figure I shows the position of 
the Rosette Nebula relative to its parent molecular cloud Mon 
OB2, as well as other known and nearby objects, marked on the 
Planck 353 GHz map (Sect. 3.1). The Rosette has been the ob¬ 
ject of extensive studies, as it has been considered an archetype 
of a triggered star formation site (Williams et al. 1994; Balog 
et al. 2007; Roman-Zuniga & Lada 2008; Schneider et al. 2010). 
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However, based on the small age spread among the young stel¬ 
lar objects in the Rosette, Ybarra et al. (2013) and Cambresy 
et al. (2013) suggest that the effect of the H ii region expansion 
on stimulating star formation is secondary relative to the origi¬ 
nal collapse of the cloud. The prominent OB association, NGC 
2244, responsible for the ionized nebula and the evacuation of its 
central part, contains more than 100 stars, of which seven are O 
and 24 are B stars (Ogura & Ishida 1981; Marschall et al. 1982; 
Park & Sung 2002; Roman-Zuniga & Lada 2008). Cox et al. 
(1990) quote a luminosity of 22 x 10^ Lq for the brightest four O 
stars and 13 B stars of the central cluster. 

The first study of the expansion of the Hit region in the 
Rosette was done by Minkowski (1949). The author also noted 
the existence of “elephant trunks” and dark globules, from the 
analysis of photographic data, in the north-western edge of the 
central cavity (Herbig 1974; Schneps et al. 1980; Carlqvist et al. 
1998, 2002). Globules are also seen in the south-eastern re¬ 
gion in the Herschel images of the Rosette (Schneider et al. 
2010), where the ionized nebula is interacting with the molec¬ 
ular cloud. The Rosette Hit region is an ionization bounded 
Stromgren sphere (Menon 1962) and expands at a velocity of 
< 15kms“* (Smith 1973; Fountain et al. 1979; Meaburn 1968). 
The mean electron temperature in the nebula, derived by Quireza 
et al. (2006), is 8500 K. 

The Rosette Hit region has a projected optical extent of 
about 1°4. Celnik (1985) estimated the size of the nebula by fit¬ 
ting several shell models to observations at 1.4 and 4.7 GHz. The 
mean radial profile of the radio emission was reproduced by a 
model of a spherical shell with inner and outer radii of 7 and 
19pc, respectively, and constant electron density of 13.5cm“^ 
(these values have been scaled from the distance of 1.42kpc 
used by Celnik to the adopted distance of 1.6 kpc). Celnik also 
found that the ionized shell shows a radial density gradient out 
to an outer radius of 28 pc, with a mean electron density of about 
5.7 cm“^. This model leads to a total mass of ionized matter of 
1.2 X 10"^ Mq, similar to the value of 1.3 x 10"^ M© obtained for 
the spherical shell of constant density. The total molecular mass 
derived from *^CO emission is 1.6 x 10^ M© (Heyer et al. 2006), 
of which 54 % lies inside the Rosette H ii region. 

The magnetic field in the Rosette nebula has been recently 
studied by Savage et al. (2013) using RM observations of back¬ 
ground galaxies. They fit the RM data with a simple stellar bub¬ 
ble model (Whiting et al. 2009) that consists of an inner low- 
density cavity of shocked stellar wind and a shell of shocked and 
photoionized ISM material. The authors assume that the ambient 
magnetic field Bq outside the bubble has a strength Bq - A jiG 
and that the field component in the shock plane is amplified by a 
factor of 4, appropriate for strong adiabatic shocks. Then by fit¬ 
ting the RM data to their model, they find that the angle between 
B() and the line of sight is 0o =72°. The results of Savage et al. 
(2013) will be further discussed in Sect. 5.2. 

Part of the northern Monoceros region has also been studied 
in polarization at radio frequencies by Xiao & Zhu (2012). These 
authors focused on the supernova remnant (SNR) Monoceros 
Loop, G205.7-1-0.1, north-west of the Rosette nebula. The SNR 
and H ii region are thought to be interacting (Davies et al. 1978; 
Smith 1973; Jaffe et al. 1997; Fountain et al. 1979; Gosachinskii 
& Khersonskii 1982). The polarization data presented by Xiao 
& Zhu (2012) do not cover the interface region between the two 
objects and show that the magnetic field in the Rosette Nebula 
is largely parallel to the Galactic plane. However, at a frequency 
of 5 GHz the observed polarization vectors are the result of ei¬ 
ther highly rotated background polarized synchrotron emission. 


considering the high RMs of about 500 rad m " measured in the 
nebula (Sect. 4.2.2, Savage et al. 2013), or foreground emission. 

3. Data 

3.1. Planck data 

Planck (Tauber et al. 2010; Planck Collaboration I 2011, 2014) 
has mapped the polarization of the sky emission in seven fre¬ 
quency channels, from 30 to 353 GHz. In this paper we use the 
data from the High Frequency Instrument (HFI; Lamarre et al. 
2010; Planck HFI Core Team 2011; Planck Collaboration VI 
2014) at 353 GHz, where the dust polarized emission is bright¬ 
est. A first description of the dust polarization sky is presented in 
a series of papers, namely Planck Collaboration Int. XIX (2014), 
Planck Collaboration Int. XX (2014), Planck Collaboration Int. 
XXI (2014), Planck Collaboration Int. XXX (2014), Planck 
Collaboration Int. XXXII (2014), and Planck Collaboration Int. 
XXXIII (2014). 

3.1.1. Intensity and polarization 

We use the full-mission temperature and polarization sky 
maps at 353 GHz from the 2015 release of Planck (Planck 
Collaboration I 2015). The map-making, calibration and correc¬ 
tion of systematic effects is described in Planck Collaboration 
VIII (2015). We smooth the maps and their corresponding co- 
variance from the initial angular resolution of 4^9 to 6'. This is a 
compromise between increasing the signal-to-noise ratio of the 
polarization data and preserving high resolution, as well as min¬ 
imizing beam depolarization effects. 

The 353 GHz intensity data are corrected for the cosmic mi¬ 
crowave background (CMB) anisotropies using the Commander 
map, presented in Planck Collaboration IX (2015), although 
the contribution to the total signal in the region under study is 
small, at most 2 %. We also subtract the Galactic zero level off¬ 
set, 0.0885 + 0.0067MJy sr ', from the intensity map and add 
the corresponding uncertainty to the intensity variance (Planck 
Collaboration VIII 2015). We do not attempt to correct for the 
polarized signal of the CMB anisotropies, as this is a negligi¬ 
ble contribution ^ 1 % at 353 GHz in the region under study 
(Planck Collaboration Int. XXX 2014). 

The linear polarization of dust emission is measured from 
the Stokes parameters Q, U, and /, delivered in the Planck data 
release. They are the result of a line-of-sight integration and are 
related as 

Q- Ip cos(2i/'), 

U - Ip sin(2i/r), 

P ^|Q^ + t/2 

- 1 -’ 

\f/ - 0.5 arctan(t/, Q), (1) 

where p is the dust polarization fraction and i/r is the polarization 
angle. The two argument function arctan(X, Y) is used to com¬ 
pute arctan(X/y), avoiding the n ambiguity. The Planck Stokes 
parameters are provided in the HEALPix (Gorski et al. 2005) 
convention, such that the angle i// = 0° is towards the north 
Galactic pole and positive towards the west (clockwise). In the 
commonly used lAU convention (Hamaker & Bregman 1996), 
the polarization angle is measured anticlockwise relative to the 
north Galactic pole. We adopt the lAU convention by taking the 
negative of the Planck U Stokes map. 
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The Stokes Q and U maps in the region under study have 
a typical signal-to-noise ratio of 8 and 2, respectively. The 
main systematic effect concerning polarization data is signal 
leakage from total intensity. We used the corrections derived 
from the global template fit, described in Planck Collaboration 
VIII (2015), which accounts for all of the leakage sources, 
namely bandpass, monopole, and calibration mismatches. The 
2015 Planck data release also includes polarization maps that 
were corrected only for dust bandpass leakage using a different 
method. We compare the two sets of maps in the region under 
study to quantify the systematic uncertainties associated with 
the leakage correction. The relative difference in the Q Stokes 
parameter is less than 10%, whereas it is typically 25% in t/ 
because the Rosette is 2° from the Galactic plane, where the 
magnetic fields are largely parallel to the plane and contribute 
more signal to Q than to U. 

Because of the presence of noise in the data, the polarization 
intensity, P, calculated directly using Eq. (1) is positively bi¬ 
ased. We debias this quantity according to the method proposed 
by Plaszczynski et al. (2014) by taking into account the diago¬ 
nal and off-diagonal terms in the covariance matrix for Q and U. 
Since we only calculate the polarization intensity towards high 
(> 10) signal-to-noise regions, namely the Mon OB2 cloud, the 
relative difference between the debiased and direct (Eq. 1) esti¬ 
mates is less than 3 %. A comparison between other debiasing 
methods as applied to the all-sky Planck data is presented in 
Planck Collaboration Int. XIX (2014). 

3.1.2. Products 


We consider 20 measurements (given in Table 3 of Savage 
et al. 2013); two of the sources are depolarized and we also ex¬ 
clude the only negative RM value detected. This measurement 
was also discarded in the analysis of Savage et al. (2013), who 
could not determine if its origin is Galactic or extragalactic. In 
any case this measurement lies outside the boundaries of the 
Rosette Nebula. Whenever two values are given for the same 
source, in the case of a double component extragalactic object, 
we take their mean and the dispersion as the uncertainty because 
we combine the RM data at a resolution of 12"8, with emission 
measure data at a lower resolution of 14!4 (Sect. 3.3). The 20 
RM measurements are positive and have a typical uncertainty of 
40 %, or below 10 % for half of the sample. 

3.3. Emission measure data 

The emission measure (EM) data are from the radio recombina¬ 
tion line (RRL) survey of Alves et al. (2014) at 1.4 GHz and at 
a resolution of 1414. These are a by-product of the Hi Parkes 
All-Sky Survey (HIPASS, Staveley-Smith et al. 1996) and their 
analysis is presented by Alves et al. (2010, 2012). The emission 
measure is defined as 



and derived from the integrated RRL emission as 

EM = 5.2xl0-(|)''’/(^)»-V. (4) 


We use the Planck “Type 3” CO map from Planck Collaboration 
XIII (2014) to qualitatively inspect the molecular gas in the 
Monoceros region. This is a composite line map, where the line 
ratios between the first three CO rotational lines are assumed to 
be constant across the whole sky. The Type 3 CO map is the 
highest signal-to-noise map extracted from Planck, with an an¬ 
gular resolution of 5.'5. 

The dust emission from Planck wavelengths to about 100 /im 
is dominated by the contribution of large dust grains (radius 
larger than 0.05/rm). In order to trace their temperature, we use 
the map derived from the all-sky model of dust emission from 
Planck Collaboration XI (2014), at a resolution of 5'. 


3.2. RM data 

We use the RM data presented by Savage et al. (2013) towards 23 
extragalactic radio sources observed through the Rosette com¬ 
plex with the Karl G. Jansky Very Large Array. ^ The RM obser¬ 
vations towards extragalactic sources are an integral of the line- 
of-sight magnetic field component By, weighted by the electron 
density iig, and given by 


RM = 0.81 



where S is the path length from the source to the observer. The 
RM is positive when B points towards us, hence when By is also 
positive. 


^ We checked for the presence of nearby pulsars in the Australia 
Telescope National Facility Pulsar Catalogue (Manchester et al. 2005). 
We did not find pulsars closer than 2° from the centre of the Rosette 
with the required information (distance, dispersion measure, RM). 


where Tg and Tl are the electron and line temperatures, respec¬ 
tively. The overall calibration uncertainty in these data is 10 %. 

3.4. IRAS data 

We use the IRIS (Improved Reprocessing of the IRAS Survey, 
Miville-Deschenes & Lagache 2005) data at 12 and 100/im. The 
IRAS map at 12/rm traces the emission from dust particles that 
are smaller than those emitting at the longer Planck wavelengths. 


4. Our observational perspective on the Rosette 

In this section we start by describing the features of the Rosette 
Nebula and its parent molecular cloud from total intensity maps 
at different frequencies. We then study the polarized emission 
arising from the dust shell that surrounds the Rosette H ii region, 
along with the radio polarization signal created by the ionized 
gas. 

4.1. Intensity 

The Rosette Nebula and its parent molecular cloud are presented 
in Eig. 2. The maps show the RRL thermal radio emission, dust 
emission at 353 GHz, integrated CO emission, dust temperature, 
and dust emission at 12/rm and 100/rm. The central cavity is 
seen both in radio and in dust emission, as a lower intensity re¬ 
gion around the position of the ionizing source NGC 2244. 

The expanding H ii region is interacting with the molecular 
cloud Mon OB2, creating photon-dominated regions at their in¬ 
terface (Cox et al. 1990; Schneider et al. 2010), which can be 
traced at 12/im. The map of Eig. 2(e) shows bright 12yum emis¬ 
sion associated with the 353 GHz contours on the eastern side of 
the nebula’s centre. Within these structures elongated along the 
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Fig. 2. The Rosette Nebula and its molecular cloud, (a) Integrated RRL emission at 1414 resolution, showing the radio morphology 
of the H II region. The green circles, with a diameter equivalent to the beam full width at half maximum (FWHM) of the RRL data, 
give the position of the 20 RM observations of Savage et al. (2013) used in this work, (b) Planck 353 GHz emission at 6' resolution. 
The black star in panels (a) and (b) indicates the position of the central star cluster NGC 2244 and the black circles correspond to 
the inner and outer radii of the Hit region, rfjj" = 7pc and = 19pc, respectively (see Sect. 4.2.2). (c) Planck CO integrated 
intensity at 515 resolution, (d) Dust temperature derived from Planck at 5' resolution, (e) Dust emission as seen by IRAS at 12//m 
and at 100/rm (f), at 4' resolution. The contours in panels (c) to (f) show the Planck 353 GHz emission; these are at every 5 % from 
15 to 30 %, at every 10 % from 30 to 50 %, and finally at 70 and 90 % of the maximum intensity of 37 MJy sr“*. 
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ionization front, there are dense molecular clumps with on-going 
star formation, as the result of the compression by the H it region. 
The dust temperature is also higher at the interface of these re¬ 
gions with the nebula, as shown in Fig. 2(d), compared to the 
lower temperature at the position of their maximum intensity. 

Cox et al. (1990) studied the Rosette complex in all four 
IRAS bands and analysed the radial distribution of the dust 
emission. These authors showed that the 12 /rm emission, which 
traces the photodissociation regions, comes from a shell sur¬ 
rounding the outer side of the ionization front, as well as from 
the molecular cloud. On the other hand, the longer wavelength 
emission at 60 and 100pm mostly arises from the H ii region and 
thus is tracing ionized and neutral gas. 

The western part of the Rosette complex has not been stud¬ 
ied as thoroughly as the main molecular cloud. It is not as 
bright in dust emission, but still visible in the Planck and IRAS 
maps (Figs. 2(b), (e), and (f)) as a fragmented shell surround¬ 
ing the central cavity. There is no significant counterpart in CO 
emission, except in the cloud located in the south-west, which 
appears elongated perpendicular to the ionization front in the 
Planck and CO maps. This cloud also has lower dust tempera¬ 
ture than the remainder of the western dust shell. 

The dust shell seen in the Planck 353 GHz map surround¬ 
ing the H II region is associated with the Rosette complex, as it 
correlates with the dust emission observed by IRAS (Cox et al. 
1990) and identified in the dust temperature map. The map of 
Fig. 2(d) shows a clear enhancement of the dust temperature in 
the region occupied by the nebula as a result of dust heating by 
the central cluster. It also illustrates a decrease in the dust tem¬ 
perature towards the central cavity because the Rosette is not 
a filled H ii region but a shell, with an inner cavity created by 
the evacuation of material due to powerful stellar winds (Smith 
1973; Fountain et al. 1979). 

The origin of the dust shell as material swept up by the ex¬ 
panding H II region will be further discussed in Sect. 5. The size 
of the shell can be estimated by fitting the Planck intensity map 
with a shell model, similar to that used by Celnik (1985) to de¬ 
rive the radius of the ionized nebula. We obtain = 18 pc and 
'"ouf “ 22 pc for the inner and outer radii of the dust shell, with a 
typical uncertainty of 2 pc. Further details are given in Appendix 
A. 

4.2. Polarization 
4.2.1. Planck data 

The maps of the Q and U Stokes parameters towards the Rosette 
Nebula and its parent molecular cloud are shown in Figs. 3(a) 
and (b), respectively. There is significant Q emission towards 
Mon OB2 relative to the background, however the polarized 
emission surrounding the Rosette H ii region does not define the 
same shell as seen in intensity. Notably, the Q and U signal in 
some regions of the dust shell is at the same level as that out¬ 
side the Rosette/Mon OB complex. Figure 4 compares the Q 
and U emission observed in the dust shell that surrounds the 
Rosette H ii region (in red) with that from an aperture that rep¬ 
resents the local background (in blue). The first distribution is 
measured within 22 pc (the outer radius of the dust shell. Sect. 
4.1) from the central cluster and the second within an aperture 
of inner and outer radii of 2° and 3°, respectively, centred at 
(/, b) - (207°, -2°) (see the circles in Fig. 3). As seen in Fig. 
4(a), the Q signal in the shell clearly deviates from that mea¬ 
sured in the background, with a higher mean value and a broader 
distribution. The distribution of U emission, shown in Fig. 4(b), 


is also broader in the dust shell, even if the mean level is similar 
to that of the background. We compare the histograms of Fig. 
4 quantitatively by means of the Kolmogorov-Smirnov statisti¬ 
cal test, which shows that the polarized signal from the shell and 
the background are significantly different. This test compares the 
cumulative distributions and thus takes into account the different 
number of points from each sample. The statistical difference 
applies when considering the dust shell as a whole, but it may 
not be valid in the case of individual structures located inside it. 
We will thus focus the analysis on the mean polarized quantities, 
without attempting to study the small-scale structure observed in 
the dust shell. This approach is also justified because our analyt¬ 
ical magnetic field model is based on a uniform density distribu¬ 
tion. 

We estimate the mean value of the /, Q, and U Stokes pa¬ 
rameters in the shell and background regions. These are listed in 
Table 1 along with the corresponding uncertainties, given by the 
standard error on the mean. The second to last column of Table 1 
lists the mean of the Stokes parameters in the dust shell corrected 
for the background emission, required for comparison with the 
model results. We note that the choice of a constant value for the 
background follows from the absence of any clear pattern in the 
polarized emission outside the Rosette/Mon OB2 that we can use 
to model the background variations. As a result, the contribution 
of the background residuals to the polarized signal of the dust 
shell needs to be taken into account. The background distribu¬ 
tions of Fig. 4 have about 8 times more points (Nback = 14 154) 
than those corresponding to the dust shell (Adust = 1 757). We 
extract Adust points from the background distribution and as¬ 
sess the variation of its mean Stokes parameters by repeating 
this exercise 10000 times. We take the standard deviation of the 
10000 mean values and add it in quadrature to the standard er¬ 
ror on the mean Stokes parameters. The final uncertainty in the 
background-subtracted values is quoted in the second to last col¬ 
umn of Table 1 . We also computed these quantities using a dif¬ 
ferent leakage correction (as explained in Sect. 3.1.1), which we 
will use in Sect. 5.3 to quantify the systematic uncertainties in 
our results. We note that the method used to compute the un¬ 
certainties does not take into account that the background fluc¬ 
tuations are spatially correlated. If we estimate the background 
in circular areas of the same radius as the Rosette dust shell, 
we find a larger uncertainty on the background-corrected val¬ 
ues. Thus, we cannot exclude that the mean Q and U values in 
Table 1 for the dust shell have a significant contribution from lo¬ 
cal background fluctuations. However, we regard the possibility 
of having a second ISM structure in the same direction as the 
Rosette as unlikely. 

We perform the same analysis in the rectangular region of 
Fig. 3 that delineates the brightest part of Mon OB2. The mean 
Stokes parameters, corrected for the background contribution 
are listed in the last column of Table 1. We use these values 
to estimate the polarization properties of the Rosette’s parent 
molecular cloud. The resulting debiased polarization fraction, 
derived according to the method introduced in Sect. 3.1.1, is 
p = (3.8+0.9)% and the polarization angle is i/i = -3° ±14°. The 
plane-of-the-sky magnetic field orientation in Mon OB2, along 
the Galactic plane, is consistent with the large-scale orientation 
seen in Fig. 3(c). 

The map of polarization fraction is shown in Fig. 3(d), where 
the lowest values of p < 3 % are seen towards the densest re¬ 
gions surrounding the Rosette Nebula. Similar values are ob¬ 
served in the less bright part of the Mon OB2 cloud at (/, b) = 
(208?5, -2?5). From the all-sky analysis of Planck 353 GHz po¬ 
larization data, Planck Collaboration Int. XIX (2014) find that 
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Fig. 3. Planck polarization maps of the Rosette region: (a) Q Stokes parameter; (b) U Stokes parameter. The dashed circles, with 
radii 2° and 3°, define the region where we estimate the background level, the full circle centred on the star cluster represents the 
outer radius of the dust shell, = 22 pc, and the rectangle encompasses the brightest part of the Mon OB2 cloud, (c) Intensity map 
at 353 GHz with the plane-of-the-sky magnetic field orientation shown by headless vectors, obtained by rotating the polarization 
angle ij/ by 90°. The vectors are plotted at every 6', from the average of Q and U within the same distance from the central pixel, (d) 
Polarization fraction p, with the same intensity contours as in Fig. 2. 


Table 1. The mean level of the Stokes /, Q, and U parameters at 353 GHz, in units of MJy sr"'. The first three columns correspond to 
the dust shell, the Mon OB2 cloud and the background regions, respectively, as defined in the text and in Fig. 3. Their uncertainties 
are given by the standard error on the mean. The last two columns list the mean Stokes parameters in the dust shell and Mon OB2 
cloud, respectively, corrected for the background contribution. The background subtraction is accounted for in the uncertainties as 
described in Sect. 4.2.1. 


Region Dust shell Mon OB2 Background Difference: Shell Difference: Mon OB2 


{I) . 9.44 ±0.08 10.6 ±0.2 3.631 ±0.006 5.80 ± 0.09 7.02 ±0.2 

(0 . 0.278 ± 0.003 0.461 ±0.005 0.1911 ± 0.0006 0.087 ± 0.004 0.270 ± 0.005 

<(/> . 0.045 ± 0.003 -0.004 ±0.005 0.0488 ± 0.0005 -0.004 ± 0.004 -0.053 ± 0.005 


higher density regions have lower polarization fraction relative 
to the maximum p^ax = 19.8%, detected at 1° resolution. This 
decrease in p can be the result of several effects, namely de¬ 
polarization in the beam or along the line of sight, variations 
in the intrinsic polarization fraction of dust grains, as well as 


changes in the magnetic field geometry. The latter was shown 
to generally explain the variations of p across the whole sky: 
Planck Collaboration Int. XIX (2014) find an anti-correlation be¬ 
tween the fluctuations in the magnetic field orientation and the 
polarization fraction. Planck Collaboration Int. XXXIII (2014) 


7 
















































Planck Collaboration: The magnetic field structure in the Rosette Nebula 




Fig. 4. Normalized distribution functions (DFs) of the (a) Stokes Q and (b) U maps measured within the background aperture (blue) 
and within the outer radius of the dust shell (red). The two regions are defined in Fig. 3 by the dashed and full circles, respectively. 


used this result to model the variations of the Stokes parameters 
across three interstellar filaments with variations of the magnetic 
field orientation, for a fixed alignment efficiency of dust grains. 
In Sect. 5 we present an analytical solution for the magnetic field 
in the Rosette H ii region and associated dust shell, for constant 
dust properties. We will use this model to explain the present 
Planck polarization observations and the radio RM data consis¬ 
tently. 


4.2.2. Radio data 


The radio emission map of Fig. 2(a) shows the Rosette Hit re¬ 
gion and its nearly circular shape. We fit the radial profile of 
the RRL emission with a shell model and find the same val¬ 
ues as Celnik (1985) for the inner and outer radii, rJjJ" = 7 pc 
and = 19 pc, respectively, with an uncertainty of 1 pc. Also 
shown in Fig. 2(a) are the positions of the RM data from Savage 
et al. (2013) (Sect. 3.2). The circles have a diameter equivalent 
to the beam FWHM of the radio survey, 14!4, in order to show 
the regions within which the electron density, n^, is estimated. 

The RRL observations provide a measure of EM for a given 
electron temperature. We use Eq. (4) with - 8500 K (Sect. 2, 
Quireza et al. 2006) to calculate EM towards the 20 RM po¬ 
sitions. The observed EM includes the contribution from the 
warm ionized gas from the background/foreground material in 
the Galaxy. We correct for this contribution by subtracting the 
average EM measured towards the RM positions located out¬ 
side the H II region. Of the 20 positions, 14 lie outside the radius 
'"ou” “ 19 pc, with an average EM of (288 + 124) cm“® pc, where 
the uncertainty corresponds to the standard error on the mean. 
The EM values measured towards the shell are in the range 
1500-5300 cm^^pc. We can thus estimate the electron density 
in the Rosette towards the remaining six positions, which are 
listed in Table 2. In the general case when the electron density 
distribution is not uniform but concentrated in discrete clumps 
of ionized gas, the filling factor / is introduced to relate the true 
path length L of a given line of sight through the nebula to the 
effective path length - fL, which is the total length occu¬ 
pied by the individual clumps. The true path length through a 
shell of inner and outer radii and respectively, is given 


Table 2. Results from the analysis of the RM and EM data: 
Column 1 gives the RM source number, as in Savage et al. 
(2013); Columns 2 and 3 list the Galactic coordinates of the 
source; Column 4 gives the linear distance from the centre of 
the Rosette; Column 5 lists the mean electron density estimated 
with Eq. (5) for f - I', Column 6 gives the RM data corrected for 
the background contribution; and the last column gives the line- 
of-sight component of the magnetic field estimated with Eq. (6) 
for f - The uncertainties in both {n^) and By are statistical, 
and thus do not include the systematic uncertainties involved, 
namely those on the distance to the Rosette, its electron temper¬ 
ature, path length (which depends on the radii), electron density 
distribution, and the calibration uncertainty of the EM data. 


Source 

1 

n 

b 

n 

[pc] 

<«e) 

[cm“^] 

RM 

[rad 

S|| 

hiG] 

16. 

205.7 

-2.1 

17.0 

9±2 

676 ± 68 

5± 1 

17. 

206.2 

-2.1 

4.6 

12 ± 1 

594 ± 30 

2.3 ±0.2 

18. 

206.8 

-2.4 

16.0 

12 ± 1 

219± 113 

1.1 ±0.6 

no ... . 

205.9 

-1.7 

16.4 

10 ± 1 

709 ± 117 

5± 1 

112 ... . 

206.6 

-1.9 

9.6 

13 ± 1 

703 ± 26 

2.1 ±0.1 

115 .... 

206.3 

-1.5 

15.9 

8± 1 

501 ± 24 

4± 1 


by (Savage et al. 2013) 



2^Hii 

^'out 


m = • 

2^Hii 

^'out 

- (f/cf 




if^ 

iff 

’ m 


where f is the linear distance between a given line of sight and 
the line of sight to the centre of the nebula, measured in the trans¬ 
verse plane through the nebula. Eor the Rosette, with —1 pc 
and = 19 pc as derived above, L(f = 0) = 24 pc is the path 
length through the centre of the shell. 

The average electron density along the line of sight is esti¬ 
mated from the EM observations using Eq. (3) as 


<«e> = 



(5) 
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Table 2 lists the average electron density at each of the six posi¬ 
tions in the Rosette for f - The derived values are between 
8 and 13 cm“^ and are consistent with those obtained by Celnik 
(1985). The mean electron density across the entire H ii region is 
12.3 cm~^, with a scatter of 4cm“^. 

Similarly to the EM data, the measured values of RM include 
the contribution from the large-scale magnetic field weighted 
by the Galactic warm ionized gas. This needs to be corrected 
for in order to study the local component of the field in the 
Rosette. The mean of the 14 measurements outside the radius 
'"ou” “ 19 pc is (132 + 20)radm“^, where the uncertainty corre¬ 
sponds to the standard error on the mean. The background mea¬ 
surements vary between 50 and 230 rad m“^ (half of these have 
an uncertainty of more than 20 %). The final six background- 
corrected RM values are listed in Table 2. 

The line-of-sight component of the magnetic field can be ob¬ 
tained by combining the EM and RM observations, using Eqs. 
(2) and (3), as 


RM/radm 

- 1 ^ mG. ( 6 ) 

0.81 V(EM/cm-6pc)(L/pc)/ 


This approximation holds if By is uniform across the nebula and 
if rif. and By are uncorrelated along the line of sight. Beck et al. 
(2003) point out that the latter assumption can lead to uncertain¬ 
ties by a factor of 2-3 in the estimated value of By. The derived 
By values in the nebula for / = 1 are listed in Table 2 and vary 
from -Hi to H-5 pG. If we consider a filling factor / = 0.1 (Herter 
et al. 1982; Kassim et al. 1989), the values increase by a factor 
of 3.2, varying between -h 3 and -h 16//G. These results are in the 
range of the By values measured in the diffuse ISM (Crutcher 
2012). As will be discussed in the next section, the hypothesis 
that the magnetic field is uniform throughout the Rosette H ii re¬ 
gion is not satisfied. The field is confined to the ionized shell and 
thus its direction must vary across the nebula. We will assess this 
aspect by means of a 2D magnetic field model and compare its 
predictions with the measured By values. 


5. The magnetic field in the Rosette 

In this section we interpret the observations by comparing them 
with the results from an analytical model of an ionized nebula 
that has expanded in a uniform and magnetized medium and 
formed a neutral shell of swept-up matter. The model is pre¬ 
sented in Sect. 5.1 and further detailed in Appendix A. We use it 
to study the RM distribution across the nebula (Sect. 5.2), as well 
as to reproduce the mean polarization of the dust shell observed 
by Planck (Sect. 5.3). 

5.1. The magnetized Stromgren shell 

The evolution of an expanding ionized nebula has been studied 
numerically, both in uniform and turbulent magnetized media 
(e.g. Krumholz et al. 2007, Arthur et al. 2011). As the H ii region 
expands, the surrounding ISM is swept up into a shell around the 
central stars. In accordance with the frozen-in condition, mag¬ 
netic field lines are dragged with the expanding gas and concen¬ 
trated in the dense shell. If the magnetic pressure is comparable 
to the thermal pressure in the H ii region, magnetic forces lead 
to departures from sphericity. Eurthermore, because the swept- 
up magnetic flux increases from the magnetic poles (along the 
direction of the initial field Bq) to the equator (90° from Bq), 
magnetic pressure in the swept-up shell tends to make the shell 
thickness similarly increase from the poles to the equator. We 


will ignore magnetic effects and assume that the H ii region ex¬ 
pands equally in all directions, creating a spherical neutral shell 
of swept-up ISM. We consider this assumption to be consistent 
with the radio observations of the Rosette Hii region, which 
show its nearly circular shape (e.g. Eig. 2(a)), despite a possi¬ 
ble elongation along the line of sight. 

Within this framework, we derive an analytical solution for 
the magnetic field in a spherical structure composed of a shell of 
swept-up gas formed around a shell of ionized gas. This config¬ 
uration is shown in Eig. A.l of Appendix A, where the details 
of the derivation are discussed. The model is an analytical de¬ 
scription of the correspondence between the initial and present 
configurations of the radial distribution of matter, which char¬ 
acterizes the expansion of the H ii region. We consider that the 
initial gas density and the magnetic field Bq are uniform and that 
they evolve as the matter expands radially. This is translated in 
the form of an expansion law (Eig. A. 2). The final magnetic field 
B (Eq. A. 8) depends on the strength of the initial field, Bq, and its 
direction, described by the polar and azimuthal angles (Oq, (po), 
with respect to the line of sight. 

We note that the Galactic magnetic field has uniform and 
random components, which are of the same order (see e.g. Beck 
2001). In the present study, Bq corresponds to the initial local 
field, which is a sample of the total (uniform plus random) mag¬ 
netic field. 

In the following sections we describe how the combination 
of RM and dust polarization data allows us to fully describe 
the magnetic field in the Rosette and its parent molecular cloud. 
These data have different resolutions; however, this is not a con¬ 
cern since the model does not attempt to reproduce any fluctua¬ 
tions on the scale of the resolution of either of the observables. 


5.2. Ionized shell: RM 

The RM of the modelled Hii region is computed using Eq. 
(2). We consider a uniform nebula with constant electron den¬ 
sity He = 12.3 cm“^, the mean value derived from the EM data 
(Sect. 4.2.2), and integrate By along the line-of-sight depth of the 
Hii shell. Since the RM is derived solely from the line-of-sight 
component of the held, it does not depend on (f>o. 

We derive the Bq and Oq values that best fit the RM data 
through a ^^-minimization procedure taking into account the 
corresponding uncertainties, as follows: 

0 (RM];nod “ EMobs) . 

^ " Zj -2-• 

A^poinK “^RM.obs 

The reduced is computed as x^ - A^/Ndof, where Ndof = 
Appoints - iVparams IS the number of degrees of freedom, the differ¬ 
ence between the number of points and the number of parame¬ 
ters in the fit. In the present case Tfparams = 2, corresponding to 
Bo and Oq. The results are shown in Eig. 5. Eigure 5 (a) presents 
the best-fit RM curve as a function of the linear distance to the 
centre of the nebula, along with the 20 RM measurements dis¬ 
cussed in Sect. 4.2.2. The blue line is the result of a fit to the 
six RM data points within the outer radius of the Rosette neb¬ 
ula, ^ = 19 pc, and hereafter referred to as the reference 

fit. The best-fit parameters are Bq = 7 pG and 6o — 25°. The 
shape of the RM curve is set by the radii of the H ii shell (thus 
by the expansion law used. Appendix A), and the scaling of the 
curve depends on a combination of the two free parameters. The 
model is unable to reproduce the scatter in the RMs observed 
close to the boundary of the H ii region, which could be due to 
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Fig. 5. Comparison between the RM observations and the pre¬ 
dictions from the model. Panel (a) shows the RM data as a func¬ 
tion of the distance to the centre of the Rosette. The two ver¬ 
tical dashed lines show the inner and outer radii of the nebula, 
rjjj" = 7pc and = 19 pc, respectively (Sect. 4.2.2). The blue 
curve shows the radial profile of the modelled RM for the best fit 
to the six RM measurements within (filled circles), the ref¬ 
erence fit. The red curve is the result of fitting the highest four 
RM data points. The open circles correspond to the RM observa¬ 
tions outside the Rosette, used to estimate the background RM, 
and are not included in the fit. Panel (b) presents the reduced 
from both fits. The stars indicate the best fit parameters Bq and 
0(), which correspond to the minimum for each fit. The refer¬ 
ence fit (blue) gives Bq = 7 jiG and 0q - 25° for;^f^ = 271 and 
Ndof = 4. The second fit (red) gives Bq — 24 fiG and 9q - 70° for 
- 76 and Ndof = 2. The contours are at 10 and 30% above 
the corresponding minimum values of;^f^. 


the possible dumpiness of the medium and/or fluctuations in the 
magnetic field direction. These effects are not accounted for in 
our model, which also assumes that the expansion is perfectly 
spherical. We did not attempt to adjust the radii of the H ii shell, 
or equivalently its expansion law. Instead, we perform a second 
fit, excluding the two lowest RM data points, at f ~ 16 pc, to 
assess the departures of the observations from the simplified as¬ 
sumptions of the present model. The result is shown by the red 
curve in Fig. 5 (a). The best-fit parameters are Bq — 24yuG and 
9q - 70°, for;t'^ = 76 with Ndof = 2, leading to aXr '^hat is about 
30 % lower than that obtained in the reference fit. 



Fig. 6. Average of the line-of-sight component of the magnetic 
field as a function of the linear distance to the centre of the neb¬ 
ula. The curve shows the result of the reference fit, for which 
Bq -1 jiG, 9o — 25°, and Boy = 6.3 //G, as indicated by the hor¬ 
izontal dashed line. The three vertical lines delineate the inner 
and outer radii of the H ii region and the outer radius of the dust 
shell: r”" = 7pc; = 19pc; and = 22pc. 

The;^'^ contours as a function of Bq and 9o for the two fits are 
shown in Fig. 5 (b). The high^^ values (given in the caption) are 
due both to the simplicity of the model and to the uncertainties in 
the RM observations, which do not reflect their true radial varia¬ 
tion across the Rosette, as they correspond to individual line-of- 
sight measurements. The contours illustrate the degeneracy be¬ 
tween the strength and the orientation of the initial field relative 
to the line of sight, as expected since RM oc By = B cos 9 (Eq. 2). 
As a consequence, all the {Bq, 9q) combinations that follow the 
minimumcontour lead to approximately the same Boy value: 
Boy = 6 ^G for the reference fit and Boy = 8 ^G for the second 
fit. Furthermore, for a given angle 9q, the two different fits give 
Bq values that differ by less than about 4 juG. 

In the rest of this section, we consider only the reference 
fit obtained when using all of the six RM measurements (blue 
curve in Fig. 5). Figure 6 shows the radial profile of the mean of 
By, (By), measured along the line of sight. There is a significant 
difference between (By) in the Hit region and in the dust shell. 
We note that the exact shape of the curve is determined by the 
adopted expansion law (Fig. A. 2), which characterizes the two 
distinct regimes in the evolution of the Rosette: the expansion of 
the ionized gas, leading to a decrease in B relative to Bq, and the 
compression of the interstellar gas in the dust shell, accompanied 
by a compression of the field lines, and hence an increase in B 
(see Fig. A. 3). The mean value of By across the projected surface 
of the ionized shell, which in 3D includes the central cavity, is 
2.6 juG. This is comparable to the mean of the six values derived 
from the RM data listed in Table 2. Therefore, the present model 
indicates that we can use the RM data to recover the mean By 
in the H ii region, which is 62 % lower than Boy in the molecular 
cloud. 

The RM observations are from Savage et al. (2013), who 
fitted a different analytical model to the data, as introduced in 
Sect. 2. The authors found a value of 9q - 72° for an assumed 
Bo = 4 juG, under the assumption that a strong adiabatic shock 
produces an enhancement of the component of the field parallel 
to the expansion front, relative to the ambient medium. In addi¬ 
tion, Savage et al. (2013) applied the shock boundary conditions 
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to the whole thickness of the shell, although these only hold in 
the thin-shell approximation. The present magnetic field solu¬ 
tion in Eq. (A. 8) naturally explains the variations of the normal 
and tangential components of the field relative to the expansion 
front, throughout the nebula. 

This analysis defines the loci of Bq and 0o values that best 
fit the RM data towards the Rosette nebula (blue curve in Fig. 
5 (b)). For all of these solutions the resulting line-of-sight field 
component in the ambient medium is Z?o|| - 6 pG, which is at the 
low end of the range of values reported by Crutcher (2012) for 
molecular clouds of similar column density as Mon OB2 (around 
3x10^^ cm“^, Planck Collaboration XI 2014). In the following 
section we show that the degeneracy between Bq and Oq can be 
alleviated by further comparing the predictions from our model 
with the Planck polarization observations towards the dust shell. 

5.3. Neutral shell: dust polarized emission 

We model the shell that surrounds the Hit region with con¬ 
stant intrinsic dust polarization fraction po and with inner and 
outer radii of - 19 pc and = 22 pc, respectively (see 
Appendix A). Since we will be comparing the ratios between 
the mean Stokes parameters, we do not specify the density or 
temperature of the gas and work with normalized quantities. The 
polarization fraction p can be written as p - po sin^ 6. We com¬ 
pute q - Po sin^ 0 cos(2i/r) and u - po sin^ 0 sin(2i/r) (Eq. 1) at 
every position in the 3D shell and integrate along the line-of- 
sight direction to obtain the normalized Stokes parameter maps. 

There are three variables involved in modelling the dust 
polarization: the angles (0o,(^o), which define the direction of 
Bq with respect to the line of sight, and the intrinsic polariza¬ 
tion fraction po. The Stokes parameters do not depend on the 
strength of the magnetic field. The polarization fraction rela¬ 
tion given above indicates that po and 0 are degenerate. We 
thus start by comparing the data with the model for a fixed 
Po value of 4%, which corresponds to the observed value in 
the Mon OB2 cloud (Sect. 4.2.1). Figure 7 illustrates how the 
predicted ratios between the mean values of the Stokes pa¬ 
rameters, {Q)I{I) and {U)I{I), vary as a function of the ini¬ 
tial magnetic field direction (0o, (po). The observed ratios, cal¬ 
culated using the values listed in the fifth column of Table 
1, are = (1.51 + 0.07 (sta.) + 0.17 (sys.)) x lO^^ and 

{U)I{I) = (-0.06+ 0.06 (sta.) ±0.41 (sys.)) x 10“^. The system¬ 
atic uncertainties correspond to the difference in the ratios when 
derived with the different leakage correction maps (Sect. 3.1.1). 
The observed ((/)/(/) ratio constrains the sky projected orien¬ 
tation of the initial magnetic field, (po ^ 0°, to within about 5°. 
This value is consistent with that measured towards the Rosette’s 
parent molecular cloud (Sect. 4.2.1) and corresponds to a mag¬ 
netic field parallel to the Galactic plane. Fixing this parameter 
allows us to study how the {Q)I{I) ratio varies as a function of 
the angle between Bq and the line of sight 0o, and the intrinsic 
polarization fraction po. This is shown in Fig. 8, for po ranging 
from 4 to 19.8%, the maximum polarization fraction observed 
across the sky (Planck Collaboration Int. XIX 2014). The com¬ 
parison between the data and the different models gives an upper 
limit on 0o of about 45° for po = 4 %. A lower limit of Oq - 20°, 
implying a field that is nearly along the line of sight, is obtained 
for the maximum intrinsic polarization fraction po = 19.8 %. 

The comparison between the present model and the Planck 
polarization observations leads to two main results. First, the in¬ 
ferred range of 0o, combined with the results from the RM study, 
restricts the range of the magnetic field strength in the Rosette’s 
parent molecular cloud to Bo - 6.5-9pG. Second, we find that 



<Q>/<I> X 10^ 

Fig. 7. Comparison between the Planck observations and the 
model. The ratios {Q)I{I) and {U)I{I) derived from the polar¬ 
ization data are given by the magenta point. The statistical and 
systematic uncertainties are shown by the solid and dashed error 
bars, respectively. The dashed grey lines show the solutions for 
(po - -5°, 0°, and 5°, with 0o varying from 10° to 80°. The solid 
grey lines indicate how the ratios change with (po, for a given 0o 
angle. The model results are calculated here for po =4%. 



<Q>/<I> X 10^ 


Fig. 8. Variation of the modelled (Q)I{I) ratio as a function of 
Po and 00, for 0o = 0° (grey lines). The observed {Q)I{I) ratio is 
shown by the vertical magenta line. 


the Mon OB2 cloud has a magnetic field structure distinct from 
that of the Perseus spiral arm; an azimuthal Galactic mean field 
with a pitch angle of -8° (Ferriere 201 1), is expected to be ori¬ 
ented at about 60° from the line of sight at the position of Mon 
OB2. Within the uniform density and polarization fraction as¬ 
sumption of this model, such values of 0o are only possible for 
a significantly low intrinsic polarization fraction po- However, 
we cannot discard depolarization effects from a turbulent field 
and/or clumpy density distribution. While several observational 
studies indicate that molecular clouds preserve the large-scale 
field orientation (see Fi et al. 2014), Planck observations ques¬ 
tion this general interpretation of almost no variation of the mag¬ 
netic fields with interstellar structures. Modelling of the dust po¬ 
larization data from Planck allows us to study the 3D geometry 
of the magnetic field in a variety of environments. The analysis 
of the magnetic field structure in nearby interstellar filaments by 
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Planck Collaboration Int. XXXIII (2014) suggests that their evo¬ 
lution is coupled to the field, which is distinct from the field of 
the clouds in which they are embedded. This result agrees with 
our findings on the Rosette/Mon OB2 complex. 

Finally, we can estimate the magnetic and thermal pres¬ 
sures in the Hit and dust shells. We do not include the possi¬ 
ble contribution from a turbulent field at scales much smaller 
than the Rosette, whose quantification is difficult. In either case 
this is expected to be small owing to the fast expansion of the 
HII region, which tends to order the magnetic field. For an ini¬ 
tial field strength Bq - 9 fiG, which corresponds to 0q = 45°, 
the modelled magnetic field has a mean value = 3.2fiG 
within the ionized nebula and = 21.4 fiG in the dust shell. 
The thermal pressure in the Hii region is ^ 2n^kT^ - 
2.9 X 10 " ergcm^^, where k is the Boltzmann constant, - 
12.3cm“^, and - 8500K (Sects. 2 and 4.2.2). The mag¬ 
netic pressure, = (B^“)^/(8^) = 0.4 x 10“'^ erg cm“^, 
is therefore smaller than the thermal pressure in the Hii re¬ 
gion, supporting our initial assumption that the Hii region is 
roughly spherical. In the dust shell, the magnetic pressure is 
^mag = 1-8 X 10“" ergcm^^, hence smaller than but of the same 
order as the ionized gas pressure. This implies that both the 
thickness and the density of the dense shell should not be con¬ 
stant (e.g. Ferriere et al. 1991). Modelling these effects is beyond 
the simplified nature of the present analysis, where we focus on 
the mean polarization properties of the Rosette, but will be taken 
into account in future work. 


6. Conclusions 

This work represents the first joint analysis and modelling of ra¬ 
dio and submillimetre polarization observations towards a mas¬ 
sive star-forming region (the Rosette Nebula) to study its 3D 
magnetic field geometry. We have developed an analytical so¬ 
lution for the magnetic field, assumed to evolve from an initially 
uniform configuration following the expansion of ionized gas 
and the consequent concentration of the surrounding ISM in a 
dense shell. The assumption of uniform density and tempera¬ 
ture distributions for both the ionized and dust shells, along with 
constant intrinsic polarization fraction of dust grains, is clearly 
an approximation. Different parts of the fragmented, swept-up 
shell presumably have distinct properties and some may even be 
pre-existing dense clouds caught by the expanding H ii region. 
Nevertheless, the model is able to reproduce the mean observed 
quantities. 

We use the Planck data at 353 GHz to trace the dust emis¬ 
sion from the shell of swept-up ISM surrounding the Rosette 
H II region. Even if the shell is clearly seen in intensity, the same 
pattern is not detected in dust polarized emission against the lo¬ 
cal background. When analysed as a whole, the polarized sig¬ 
nal from the dust shell is significantly distinct from that of the 
background and can be reproduced by the current model. The 
correspondence between the model and the Planck observations 
constrains the direction of the magnetic field in the Rosette’s par¬ 
ent molecular cloud Mon OB2 to an angle in the plane of the sky 
(po ^ 0° (roughly parallel to the Galactic plane) and an angle to 
the line of sight 6 < 45°. This result is crucial to removing the 
degeneracy between Oq and Bq inherent in the RM modelling. 
We thus find that Bq is about 6.5-9/iG in Mon OB2. 

The present magnetic field model provides a satisfactory fit 
to the observed RM distribution as a function of the distance 
from the centre of the Rosette H ii region. More data are needed 
to better understand the abrupt variations of RM close to the 


outer radius of the nebula. The RM modelling suggests a signif¬ 
icant increase in the line-of-sight magnetic field from the H ii re¬ 
gion to the dense shell, where By reaches nearly 4 times Boy- 
The combination of RM and dust polarization data in this 
work is essential to constrain both the direction and the strength 
of the field in the Rosette region. The model presented here can 
be directly applied to other similar objects for which the expan¬ 
sion law can be derived. 
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Appendix A: Magnetic fieid modei 

In this appendix we present the derivation of the analytical for¬ 
mula that describes the magnetic field structure in a spherical 
shell, following the expansion of an ionized nebula in a uniform 
medium with density no and magnetic field Bq. 

We start by deriving the expansion law, which will define 
how the initial uniform magnetic field is modified. Once the star 
cluster is formed, it ionizes the surrounding gas, which becomes 
overpressured and starts expanding at a velocity close to the 
ionized gas sound speed. The expansion of the ionized gas, in 
turn, creates a cavity and sweeps up the surrounding ISM into a 
thin and dense shell. The resulting structure, as observed at the 
present time, is composed of a cavity of radius surrounded 
by a thick shell of ionized gas extending from to it¬ 
self surrounded by a thin dust shell extending from to 
(see Fig. A.l). For the Rosette, rfjJ" = 7pc and = 19 pc, as 
derived from the radial distribution of the radio emission (Sect. 
4.2.2); = 18 pc and = 22 pc, measured from the Planck 

353 GHz latitude cut through the centre of the shell (Sect. 4.1). 
Thus, the inner radius of the dust shell is slightly smaller than 
the outer radius of the H ii region. This is not surprising, as the 
Rosette is an ionization bounded Stromgren sphere and thus the 
two shells are expected to overlap at the boundary, where the 


ionized and neutral gases are mixed. However, for the sake of 
simplicity, we take = 19pc. 

We use a parameter grid of 1 pc resolution for all the radii. 
Therefore, we adopt an uncertainty of 1 pc in and We 

ascribe a larger uncertainty of 2 pc to and ; this reflects 

the variation of the fitted radii when considering the longitude 
or radial profile of the dust emission, which are affected by the 
presence of the Mon OB2 cloud. 

We now denote by ro the initial radius of a particle currently 
at radius r. In our simplified model, the initial radius of a par¬ 
ticle currently at rjjj" is simply ro - 0, while the initial radii 
of particles currently at and can be denoted by " 
and respectively. We emphasize that " and are just 
two working quantities, which do not correspond to any physical 
boundaries. The value of can be obtained by noting that a 
particle currently at has just been reached by the expanding 
shell and has not yet moved from its initial radius so that 
^dust _ ^dust _ 22 pc. The value of “ can be inferred from the 
conservation of mass. In the initial state 

j^(ro''fno = Mhu, (A.l) 

p ((r^')^ - (r«'')3) no = Mou.t, (A.2) 


where Mhh is the mass inside the sphere with radius “ and 
Mdust is the mass inside the shell with inner and outer radii " 
and respectively, as shown in the left panel of Fig. A.l. 
From the previous equations we can write 


„dust 


^dus 

4 ^ 


-I- 1 


- 1/3 


(A.3) 


We calculate the mass of ionized gas based on the size of the 
nebula and on the mean electron density of 12.3 cm“^, derived 
from the radio data (Sect. 4.2.2). Taking into account the con¬ 
tribution from ionized helium (Celnik 1985), we obtain Mhh = 
1.2 X 10^ Mq. For the mass of the dust shell we use the results 
of Heyer et al. (2006) (Sect. 2), Mdust = 8.6 x 10"^ Mq. As a re¬ 
sult, Lg “ = 11 pc. With the three known radii in the initial state, 
(0, and their corresponding values in the present state, 

Lqu”, Lpjjf), we can derive an expansion law as shown in 
Fig. A.2. With only three data points, we consider the simplest 
description of the expansion by writing r as a piecewise linear 
function of lq. This is in any case sufficient to describe the two 
clear regimes seen in Fig. A.2; the expansion of the Hii region 
(slope larger than 1) and the compression of the ISM within the 
outer shell (slope smaller than 1). The expansion law is hence 
given by r = aro + [5, where a = 1.09, yS = 7 pc for the H ii re¬ 
gion (ro < Lq “) and a - 0.27, yS = 16 pc for the dust shell 
(rH'' < ro < r^-O. 

Now consider a cartesian coordinate system with z-axis 
along the line of sight and x- and y-axes in the plane of the sky, 
along the trace of the Galactic plane and along the rotation axis, 
respectively. The initial uniform magnetic field can be written in 
these cartesian coordinates as 


with 


Bo — Box^x T BoyCy -H Boz^z^ (A.4) 


Bqx - Bo sin 6o cos (po 
Boy - Bo sin 6o sin (po 

Boz - Bo cos 00, (A.5) 
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Fig.A.l. Sketch of the adopted spherical configuration. Left the initial state, with uniform density and magnetic field; there are no 
separate ionized and dust shells. Right present-day state with a cavity inside an ionized shell between and and a dust 
shell between and . The present-day radii and correspond to initial radii 0, and respectively. 



Fig.A.l. Radial expansion law derived from the three known 
radii in the initial, ro, and final, r, states (red filled circles). The 
black curve is described by a piecewise linear function of the 
form r - aro + y6, with a = 1.09, J3 - 7 pc for ro < r^" and 
a = 0.27, j8 = 16pc for r«" < ro < r;j“‘. 

where ffg and fo are the polar and azimuthal angles, respectively, 
i.e. 00 is the angle between Bq and the line of sight, while fo 
gives the plane-of-the-sky direction of Bo with respect to the 
Galactic plane (fo - 0° for Bo parallel to the Galactic plane). 
The vector potential associated with Bq is 

Ao = Boyzejc + BozXCy + Boxye^, (A.6) 

which satisfies the condition Bq -77 x Ao. We use the frozen-in 
approximation and assume that the magnetic field evolves from 
the initially uniform configuration following the radial expan¬ 
sion of the gas. The vector potential in the final state is given by 
(see equations 4 to 10 in Parker 1970) 

A(r) = (Vr„) ■ Ao(ro), (A.7) 

and the resulting magnetic field works out to be 

j Bor^r +-^^^{Bogeg + Bo^ef). (A.8) 

The previous equation, written in spherical coordinates, clearly 
shows the change in both the normal (radial) and tangential com¬ 
ponents of the magnetic field, relative to the expansion front. 


We create a 3D cartesian grid of 181^ = 5 929 741 voxels, 
each equivalent to 0.25 pc (0'5 at the distance of the Rosette), 
with the bubble-shell structure located at the origin. The resolu¬ 
tion of the model is finer than that of the observations, which is 
needed to have the required sampling to compute the integrals 
along the line of sight. We use Eq. (A.8) along with the expan¬ 
sion law of Fig. A.2 to calculate the magnetic field strength in 
every pixel of the grid. Figure A.3 shows how the field strength 
in the shell, B, varies relative to the initial field strength, Bq. The 
map corresponds to a vertical cut through the centre of the shell 
for an initial field with (0o, fo) - (90°, 0°), therefore on the plane 
of the sky and along the Galactic plane. Figure A.3 illustrates 
that the largest compression of the field lines occurs towards the 
equator of the shell, or in the direction perpendicular to the ini¬ 
tial field Bo, where the ratio BjBo is seen to increase from the 
centre to the outer radius of the dust shell. The change in ex¬ 
pansion law at the boundary between the Hit and dust shells, 
r - 19 pc, results in a discontinuity in the tangential component 
of the magnetic field. On the other hand, close to the poles of 
the shell, or along Bo, the field lines are little disturbed, with the 
ratio BjBo continuously increasing from the centre to Z?/Z?o = 1 
at the boundary of the dust shell. Owing to the axial symmetry 
of the magnetic field model, the map of Fig. A.3 is reproduced 
in every plane about the direction of the initial field Bq. 
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Fig.A.3. Ratio of the field strength B/Bq for a vertical slice 
through the centre of the bubble-shell structure {xy plane). The 
initial field Bq is on the plane of the sky (6q - 90°) and along 
X {(po - 0°). The three black circles delineate the radii of the 
Rosette Hit and dust shells: = 7 pc; = 19pc; and 

Cf = 22pc. 
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